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Abstract 


A structural and micro-pore analysis of a series of heat treated electrolytic manganese dioxide (EMD) samples has been conducted. In terms 
of crystal structure, the original EMD with y-MnO;, structure (orthorhombic unit cell) was found to progressively convert to B-MnO, (tetragonal 
unit cell) at elevated temperatures. The structural transition was kinetically limited, with the higher temperatures leading to a greater degree of 
transformation. The orthorhombic y-MnO, unit cell was found to contract along the a and b axes, while along the c axis an expansion was observed 
only at the highest heat treatment temperatures. These changes occur as a result of manganese ion diffusion leading to the formation of a denser, 
more defect free material. The porosity of these heat treated EMD samples was also examined by Nz gas adsorption combined with various 
interpretive techniques such as the Kelvin equation, MP method, Dubinin—Radushkevich method, Dubinin—Astakhov method and a more modern 
density functional theory based approach. Despite shortcomings associated with certain techniques, all clearly indicated that the EMD micro-pore 
volume decreased and the meso- and macro-pore volume increased as the heat treatment temperature was increased. This was justified as a result 
of manganese ion movement during the structural rearrangement causing the small pores to be progressively sintered shut, while the larger pores 


were formed as a result of stress-induced cracking in the denser final product. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 
1.1. Background 


Manganese dioxide has for over 50 years been used as the 
cathode active material in primary batteries. This widespread 
popularity is due to a unique combination of its physico- 
chemical and electrochemical properties, including its relatively 
high density and purity, its ability to adequately maintain up to 
and including intermediate discharge rates for prolonged periods 
of time (particularly when used in concentrated alkaline elec- 
trolytes), as well as the fact that it is environmentally benign 
when disposed. Add to that the inexpensive nature of the raw 
materials necessary to prepare manganese dioxide, leads to it 
being a very difficult material to supercede. Despite many years 
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of usage and research into its behaviour there are still many 
properties of manganese dioxide that are unclear and need clari- 
fication. One of these is the role of manganese dioxide porosity, 
which is the subject of this work. 


1.2. Manganese dioxide structures and morphology 


As a compound, manganese dioxide has a wide array of 
different crystal structures, so much so that over 30 different 
species have been classified [1]. The importance of structure is 
that it influences dramatically the electrochemical behaviour of 
the material. It is therefore essential that we describe in some 
detail the structures we will encounter in this work. 

Manganese in oxidation states II, II and IV prefers almost 
exclusively an octahedral coordination, albeit somewhat dis- 
torted. This is particularly true for Mn(II) containing species, 
since it is a Jahn—Teller distorted ion. Therefore, the basic build- 
ing block for the majority of manganese oxides and hydroxides 
is the [Mn(O?~,OH~,OH2)6]** octahedral unit. Various cor- 
ner and/or edge sharing assemblages of these octahedral units 
give rise to different structures [2]. An alternative approach 
to interpreting manganese oxide structures is to visualize a 
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hexagonal close-packed array of oxide anions, into which man- 
ganese cations are placed selectively into the octahedral sites. 

The simplest manganese dioxide structure is that of B-MnO» 
(pyrolusite) which consists of corner sharing octahedra in the 
a-b plane, and edge sharing octahedra in the c direction. The 
result is a [1 x 1] tunnel structure progressing in the c direction, 
as shown in Fig. 1(a). For this tetragonal structure the unit cell 
parameters are ag = 0.439 nm and co = 0.287 nm [1]. In a similar 
way, the structure of ramsdellite can be visualized, this time 
though with both edge and corner sharing in the a—b plane to 
form a [2 x 1] tunnel headed into the c direction. Ramsdellite 
has an orthorhombic unit cell with ag = 0.446 nm, bp = 0.932 nm 
and cp =0.285 nm [1]. A schematic of its structure is shown in 
Fig. 1(b). 

The preferred phase for use in the alkaline manganese dioxide 
cathode is y-MnOz. The basis for describing the y-structure is as 
a random microscopic intergrowth between both the pyrolusite 
and ramsdellite structures of manganese dioxide [3], as shown 
in Fig. 1(c). Other features associated with the y-MnOz struc- 
ture include cation vacancies, lower valent manganese cations 
(Mno), and structural water present as protons associated 
with oxide anions so as to compensate for the charge deficiency 


Fig. 1. Schematic representations of: (a) the pyrolusite (8-MnOz), (b) the rams- 
dellite, and (c) the y-MnOz structures. 


incurred by cation vacancies and Mn(III) ions [4]. As a result, 
the chemical composition of y-MnOz has been described as 


(Mn**),—.—~y-(f)x-(Mn?*)y-(077)2—4x-y-(OH arty (1) 


where x and y represent the mole fraction of cation vacancies 
(O) and MncII) ions, respectively. Typical ranges for x and y 
in an unreduced y-MnO7 sample are 0.06—0.08 and 0.04-0.12, 
respectively. Another structural component of -y-MnO> that has 
been proposed is microtwinning in the (0 2 1) and (06 1) growth 
planes. Structural modeling including microtwinning has been 
used to explain the appearance of certain features in the X- 
ray diffraction pattern of y-MnOz; however, there is conflicting 
physical evidence as to its actual existence [5,6]. 

Clearly, based on the currently accepted structural and chem- 
ical (Eq. (1)) definition, there can be a multitude of materials 
loosely classified as y-MnObo, even though their structural, and to 
some extent chemical, composition can vary considerably. There 
are many synthetic routes leading to y-MnOz, each with their 
own list of experimental variables that can alter its structural and 
chemical composition as defined above. However, the y-MnO2 
used as the cathode-active material in alkaline cells is most often 
prepared using electrolytic means (electrolytic manganese diox- 
ide, EMD) [7], which in itself imparts a unique morphology to 
the material. The original high resolution TEM work on EMD 
was conducted by Turner and Buseck [8], and it showed defini- 
tive evidence in support of De Wolff’s microscopic intergrowth 
model for y-MnO> [3]. More recent TEM work on EMD by 
Heuer et al. [9] has shown the material to consist of needle-like 
crystallites ~10nm x 10nm x 100nm in size. These needles 
were then found to be assembled in a somewhat ordered fashion 
into grains ~200 nm in diameter. These grains were then found 
to be arranged randomly to form the bulk structure. Also appar- 
ent between the crystallites and grains were pores, into which 
electrolyte is expected to permeate during electrochemical dis- 
charge. It is these pores, their size and distribution, that we intend 
to characterize in this work. 


1.3. The importance of porosity 


Before discussing manganese dioxide porosity it is impor- 
tant to emphasize why it is so important. As mentioned above, 
the structural variety of manganese dioxide preferred for use in 
alkaline Zn/MnOy) cells is y-MnO . As a cathode active material, 
y-MnO> undergoes electrochemical reduction according to 


MnOOH, + AH20 + Ae — MnOOH,+4 + OH™ (2) 


where r is the mole fraction of hydrogen (H*/e~ pairs) present 
in the starting material (7 ~ 0.1 in an unreduced sample), and A 
indicates the mole fraction of hydrogen inserted during reduction 
[10]. This is a homogeneous reduction process for the major- 
ity of the discharge in concentrated alkaline electrolytes; i.e., 
down to r+ A ~ 0.8 [11]. The reduction itself involves electron 
injection from the external circuit (via the conductive graphite 
chain) into the manganese dioxide structure to reduce Mn(IV) 
to Mn(II]), together with proton insertion at the manganese 
dioxide-electrolyte interface as a result of water decomposition. 
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Fig. 2. Schematic of the y-MnO% reduction mechanism. 


Once inserted, the Ht/e~ pairs diffuse under an activity gradi- 
ent to achieve a uniform distribution throughout the structure 
(Fig. 2). Given that proton transport is widely believed to be 
the rate limiting property in the discharge of manganese dioxide 
[12], the facility with which protons are inserted into the struc- 
ture is a key electrochemical property. Inherent with this is of 
course the size of the manganese dioxide-—electrolyte interface, 
which is for the most part associated with pores. 


1.4. This work 


Our approach in this study has been to examine the changes 
in porosity that occur when a sample of electrolytic manganese 
dioxide (EMD) is subjected to a variety of heat treatment tem- 
peratures. Our intention is to begin to understand how physical 
changes in the structure of the material influence its morphology. 


2. Experimental 


2.1. Material synthesis and characterization 


The starting material used in this work was a standard com- 
mercial EMD provided by Delta EMD Australia Pty Limited. 
It was prepared by anodic electrolysis of a hot (~98 °C), acidic 
(H2SO,4) solution of MnSO, onto a titanium substrate. After 
deposition the EMD was mechanically removed from the anode 
substrate, milled into a —105 zm powder, washed, neutralized 
and then dried before being ready for use. 

Sub-samples of this starting material were then taken and 
subjected to heat treatment in air at temperatures from 100 up 
to 400°C (in 100 °C increments). The duration of the heat treat- 
ment was 24h. 

Structural information was obtained on the materials using 
X-ray diffraction (XRD). A Philips PW1710 diffractometer 
equipped with a Cu Ka source was used. Patterns were recorded 
over the 26 range 10-80° with a step size of 0.05° and a count 
time of 2.5 s. 

Morphological information on each sample was obtained 
through the use of a Philips XL30 scanning electron microscope 
(SEM). 


2.2. Porosity determination 


Sample porosity was determined from gas adsorption/ 
desorption isotherms collected using a Micromeritics ASAP 


2010 Surface Area and Porosity system, adapted with a turbo- 
molecular pump to examine low pressure adsorption. Prior to 
analysis, samples were out-gassed under vacuum at 110°C for 
2h to remove surface water which would otherwise interfere 
with the measurement. N2 adsorption and desorption isotherms 
were then collected at 77 K. 


3. Results and discussion 
3.1. Crystal structure and morphology 


X-ray diffraction is the most common tool used to character- 
ize the crystal structure of manganese dioxide. However, despite 
its widespread use, various forms of manganese dioxide tend to 
be finely divided and disordered meaning that their XRD patterns 
often contain broad peaks which overall makes the technique of 
powder XRD less useful. Nevertheless, the structural modeling 
work conducted by Chabre and Pannetier [5] does allow us to 
shed some light on the structure of EMD from its XRD pattern. 

XRD patterns of the samples prepared in this work are shown 
in Fig. 3, with the pattern for EMD labeled with the correspond- 
ing Miller indices assuming an orthorhombic unit cell [5]. The 
most obvious changes occurring as a result of heat treatment 
are (i) merging of the (1 1 0) (~22°26) and the (1 3 0) (~36°260) 
lines in the EMD (y-MnOz) starting material to form a single 
peak at ~29°26 for the 400°C heat treated sample; (ii) clearer 
separation of the (22 1) and (240) lines at ~56°20 by 400°C; 


Normalized Intensity 


EMD 


20 


Fig. 3. XRD patterns of heat-treated EMD. 
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and (iii) disappearance of the line at ~68°26, with the corre- 
sponding emergence of two lines at ~65°26 and ~73°26. These 
observations concerning the XRD pattern are consistent with 
previous works [13] in the sense that they indicate the progres- 
sive conversion of y-MnOz2, with a predominantly ramsdellite 
structural composition, into pyrolusite or B-MnOz. To quantify 
this change, the approach developed by Chabre and Pannetier 
[5] for XRD pattern analysis has been adopted to determine the 
fraction of pyrolusite in the y-MnOp structure (P,), as shown 
in Fig. 4. What is most significant is that even at relatively low 
temperatures, and after only 24h at temperature, the y-MnO2 
structure has already begun to convert to pyrolusite. 

This can be reinforced further by considering the changes 
in unit cell parameters (extracted from the XRD patterns) as a 
function of temperature. For the EMD sample used in this work 
its unit cell parameters were determined to be ag = 0.4475 nm, 
bo = 0.9676 nm and co = 0.2825 nm, which is slightly expanded 
in the a-b plane, but contracted along the c axis compared 
to ramsdellite. The most convenient way to demonstrate how 
the structure changes with heat treatment was to determine the 
difference between the unit cell parameters of the heat treated 
samples and the starting EMD, as shown in Fig. 5. As shown, 
both the ao and bo parameters decrease progressively, essen- 
tially from the outset of heat treatment. However, co remains 
relatively constant throughout most of the temperature range 
studied, increasing substantially only at 400°C. If we consider 
the respective y- and B-MnOy} structures, as shown in Fig. 1, the 
main differences are apparent in the a—b crystallographic plane. 
The fact that the unit cell contracts along the a and b axes sug- 
gests that ion (Mn**) movement in these directions predominates 
during heat treatment, so as to ultimately form a denser, more 
defect free structure. Along the c axis, the structural expansion 
at 400 °C may be due to an excess of edge sharing octahedra in 
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Fig. 4. Pyrolusite fraction (P,) within the heat treated EMD samples. 
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Fig. 5. Changes in unit cell parameters as a function of heat treatment temper- 
ature. 


a uniform array, without the presence of defects such as cation 
vacancies to compensate or provide a buffer for the close prox- 
imity of the Mn** ions. 

SEM images of the 100 and 400 °C heat treated EMD samples 
are shown in Fig. 6. These images are of the surface of a much 
larger particle, and demonstrate that both samples are composed 
of agglomerates of small grains less than ~ 100 nm in size, which 
is still much larger than the average crystallite (primary particle) 
size calculated from the Scherrer equation [14] using the XRD 
data in Fig. 3; i.e., ~14 nm, suggesting that these grains are also 
agglomerates. Itis also apparent that there is very little difference 
in morphology between these two samples, indicating that heat 
treatment has little effect on morphology in this size range. 


3.2. Gas adsorption/desorption isotherms 


An example of a typical gas adsorption and desorption 
isotherm collected in this work is shown in Fig. 7. Both the 
adsorption and desorption isotherms were as expected, particu- 
larly the desorption isotherm which demonstrated a substantial 
hysteresis loop in the P/Po range above 0.4. This is believed due 
to the slower than expected evaporation from the surface as a 
result of surface tension phenomena associated with the lique- 
fied adsorbate [15]. Furthermore, at relative pressures greater 
than 0.9, the rapid increase in the volume of adsorbed gas 
is due to condensation of the adsorbate in the larger “pores” 
formed between particles. This has been shown to be the case 
as a result of studies on different particle size distributions of 
the same material [15]. At very low relative pressures below 
P/Po =0.1, the adsorbate shows a high affinity for the solid 
surface, approaching zero only at relative pressures of ~10~°. 
Needless to say, gas adsorption in pores contributes substantially 
to the adsorption behaviour at these relative pressures. 
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Fig. 6. SEM images of EMD heat treated at: (a) 100°C and (b) 400°C. 


The first analysis of the data that can be conducted is deter- 
mination of the surface area using the linearized BET isotherm 
[16]; i.e. 

x o 1 
Val =x)  VmC 


(C — 1)x 
VinC 


(3) 


where x= P/Po, Va the volume of gas adsorbed at x (cm? g7!), 
Vm the volume of gas that would occupy a monolayer (cm? g7!), 
and C a constant for the gas/solid pair. From a plot of x/V,(1 — x) 
versus x, Vm can be determined from which the surface area 
(S; m? g7!) can be calculated using the expression S=4.35Vm, 
in which case the factor of 4.35 takes into account the cross- 
sectional area of the adsorbing N2 molecule (16.2 A?) and vari- 
ous other conversion factors. The BET isotherm can only really 
be considered valid over the partial pressure range from 0.05 
to 0.30 because at low partial pressures the isotherm under- 
predicts actual adsorption due to higher than expected gas 
adsorption energies on the solid as a result of surface inho- 
mogeneity, while at relatively high partial pressures the BET 
isotherm over-predicts adsorption due to limitations associated 
with liquefication of the adsorbate [15]. BET surface areas of 
the heat-treated EMD samples, determined using Eq. (3), are 
shown in Fig. 8 as a function of the heat-treatment temperature. 
Firstly, the surface area of these samples are relatively high, 
much greater than the geometric surface area predicted from 
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Fig. 7. Gas adsorption and desorption isotherms for the starting EMD used in 
this work. (a) Linear and (b) logarithmic P/Po data. 


the particle size distribution. This of course indicates that the 
materials are quite porous. Secondly, note that the calculated 
surface area decreases progressively as the heat-treatment tem- 
perature is increased, particularly at the highest temperature, in 
which case the surface area decreases by ~50%. Given the struc- 
tural changes that are occurring upon heat-treatment, with the 
corresponding diffusion of Mn** ions through the oxygen anion 
framework, the decrease in surface area suggests that EMD pores 
are being sintered shut, essentially further crystallizing the mate- 
rial structure. This will of course be confirmed by examining the 
changes in pore size distribution with heat treatment. 

There are of course alternative isotherms, other than BET 
theory, that can be used to describe gas adsorption on a solid 
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Fig. 8. BET surface areas of samples prepared in this work. 


substrate. One such approach is the Frenkel, Halsey and Hill 
(FHH) isotherm [17]: 


b 
i fa : 
m n(Po/P) 


where a and b are empirical constants determined by fitting the 
isotherm to the experimental data. Another is the Harkins and 
Jura (HJ) isotherm [18]: 


Va d! 0.5 
Vin b' — log(P/ Po) 


where a’ and b’ are again empirical constants. What is appar- 
ent about these two isotherms is that Vm cannot be determined 
directly, and hence the surface area cannot be estimated from 
these isotherms alone. However, if we were to employ the val- 
ues of Vm obtained from the BET analysis, then estimates of 
the empirical constants a and b (or a’ and b’) can be made by 
fitting either Eq. (4) or (5) to the experimentally obtained adsorp- 
tion data (Table 1). The partial pressure range over which the 
isotherms were fit was 0.05—0.85, and as can be seen in Table 1 
from the 7? values, the fit over this range was excellent in both 
cases. Outside of this partial pressure range significant devi- 
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Fig. 9. Example of a t-plot using the FHH and HJ models (Eqs. (6) and (7), 
respectively). 


ation between the experimental and predicted isotherms was 
observed. 

The significance of the FHH and HJ isotherms is that they 
allow us to determine the statistical thickness of the adsorbate (£), 
which in turn can be used to estimate the pore size distribution 
of the solid. Knowing V,/Vm, which is the relative volume of gas 
adsorbed compared to the monolayer volume, the layer thickness 
can be determined using: 


b 
a o 
d 0.5 
HJ t values: t= 354 EFES (A) (7) 
b! — log(P/ Po) 


where 3.54 A corresponds to the thickness of a N3 layer. An 
example of such a t-plot is shown in Fig. 9, which uses the 
absorption data for the starting EMD. Of course the data in Fig. 9 
does not have much significance below a t value of 3.54 A since 
this is the thickness of a N2 monolayer. However, what it does 
indicate is that there is only partial coverage of the surface with 
N2 molecules. Above 3.54 A we have greater than a monolayer 
coverage, with the t value indicating the average N2 thickness. 


Table 1 
Empirical parameters for the FHH and HJ isotherms 
Sample Vm (cm? g7! at STP) FHH HJ 

a b r da b' r? 
EMD 7.70 4.011 0.267 0.982 1.513 0.242 1.000 
100°C 7.06 3.819 0.295 0.971 1.448 0.185 1.000 
200°C 6.74 3.079 0.362 0.978 1.239 0.093 0.997 
300°C 6.58 2.717 0.428 0.991 1.133 0.037 0.998 
400°C 4.05 2.196 0.548 0.993 0.968 —0.018 0.998 
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3.3. Pore size distributions 


Many different methods for assessing the porosity of solid 
materials from gas adsorption and desorption data have been 
proposed. Methods specific for meso- and macro-porous mate- 
rials, as well as micro-porous materials have been reported and 
used extensively. In general, these methods of analysis are quite 
empirical in nature, and in some instances have flawed assump- 
tions, but nevertheless they still find widespread use despite these 
shortcomings. 


3.3.1. Using the Kelvin equation 

For porous materials the Kelvin equation has been used exten- 
sively despite being relevant for only meso- and macro-porous 
materials [15,19]. The premise underlying this method is that 
the vapour pressure of the adsorbate in equilibrium with con- 
densed adsorbate in the pores is less than Po as on a flat surface. 
On changing the external pressure, gas will adsorb or desorb 
from the available free surface, and when a suitable pressure is 
reached the centre of the pore will either fill from an adsorbed 
layer on the pore walls, or empty leaving behind an adsorbed 
layer. For cylindrical pores, the radius of pores is given by the 
sum of the statistical thickness t, as determined using thickness 
equations such as in Eqs. (6) and (7), and the Kelvin radius; i.e. 


e(2)- 


where y and Vj, are the surface tension and molar volume of 
the condensed liquid adsorbate, 6 the contact angle between the 
solid and condensed phase (for Nz as the adsorbate, @=0° and 
so cos = 1), rm is the radius of curvature of the meniscus, while 
R and T have their usual significance. 

Using such an approach, the calculated pore size distribution 
for our heat treated series of EMDs is shown in Fig. 10. Using this 
method it is not straightforward to extract the behaviour of the 
micro-pores (<20 A) as the initial EMD is heat treated. In some 
ways this was expected given that the use of the Kelvin equation 
is most relevant for meso-pores (20-500 A) and macro-pores 
(>500 A). Above 20 A the resolution of the pore size distribution 
becomes sufficient such that we can extract meaningful informa- 
tion. What Fig. 10 suggests is that with heat treatment of EMD 
the pore volume (Vp) increases at larger pore sizes as the heat 
treatment temperature increases. As suggested previously, the 
reason for this is as the heat treatment temperature increases, 
the mobility of manganese ions in the structure increases to the 
point where the smaller pores within the structure are being 
sintered shut as a result of ion movement. However, with ion 
movement comes a structural rearrangement from y-MnOz to 
the more dense B-phase. This increase in density may causes a 
stress build up in the heat treated materials, which in turn induces 
the formation of more larger cracks (larger pores) within the 
material. 

One of the objectives of this work was to compare various 
methods of pore size analysis to confirm the trends we observe 
using any particular method. To do this correctly we must under- 
stand the various pitfalls associated with each method. The 
Kelvin method has been found to break down when P/Pg reaches 


2yVi cos 6 
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Fig. 10. Pore size distribution determined using the Kelvin equation combined 
with HJ thickness (f) data. 


about 0.4 because the liquid nitrogen cannot support the high 
pressure differential of over 0.5 atm [20]. As aresult, many mate- 
rials appear to have a lot of pores at about 40 A, when in fact it is 
a breakdown of the assumption of liquid behaviour that is being 
measured and actual pore sizes may be significantly smaller. For 
the data presented in Fig. 10, the peak in pore volume occurs 
around 40 A bringing into question the validity of the calculated 
distribution. However, a redeeming feature of our data is the fact 
that an adsorption isotherm was used, which in effect avoids the 
hysteresis loop apparent in the desorption isotherm that gives 
rise to the overestimation of 40 A pores. 

While use of the Kelvin equation focuses primarily on larger 
pore sizes, as is clear from the results in Fig. 10, alternative meth- 
ods have been developed to examine the smaller micro-pores. 
As we shall see these methods are quite empirical in nature, and 
as such the quality of results can be quite limited. This is par- 
ticularly true for micro-porous materials where pore sizes are 
comparable to adsorbate dimensions. 


3.3.2. MP method 

The MP method of pore analysis [21] also makes use of the 
adsorbate thickness data described previously. From a plot of 
volume of gas adsorbed (V,) versus the statistical thickness (t) 
the pore size distribution can be extracted by determining the 
equation of the tangent at each t value. The value of the tangent 
slope corresponds to the total micro-pore area in unfilled pores 
up to that t value. The pore volume (Vp) at this point can be 
determined from: 


(Sn = Snti(tn F tn+1) K 


Vp = z 


(9) 


where S; is the slope of the tangent at thickness t;, and K is 
a constant converting gas volume at STP to liquid volume, 
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Fig. 11. Pore size distribution calculated using the MP method and FHH thick- 
ness data. 


and Å to cm. The data in Fig. 11 show the results of using 
such an approach (combined with the FHH thickness data) 
to determine the pore size distribution. Similar trends can be 
observed with the use of HJ thickness data. In many ways 
the data in Fig. 11 confirms our model of structural change. 
Namely, as the heat treatment temperature increases the pore 
volume associated with very small pores decreases, as seen in 
the data at ~5 A. Furthermore, the pore volume associated with 
slightly larger pores (~7 Å in the initial EMD, up to ~15 Å for 
the sample with the highest heat treatment temperature) also 
increases with heat treatment temperature. Compared to using 
the Kelvin method, these trends in pore size distribution are sim- 
ilar but cover different pore size ranges, and so are not directly 
comparable. Nevertheless, the similarities do give us some 
confidence. 


3.3.3. Dubinin—Radushkevich (DR) and Dubinin—Astakhov 
(DA) methods 

While the concept of multi-layer gas adsorption on flat sur- 
faces and the walls of large pores is valid, it is not neces- 
sarily appropriate when the pores themselves are comparable 
in size to adsorbate molecular dimensions. This led Dubinin 
and Radushkevich [19,22] to develop the following empirical 
expression for the filling of pores in microporous carbons; i.e. 


A 2 
rapel (=) (10) 


where V is the volume adsorbed at P/Pp and temperature T, Vo 
the limiting micro-pore volume, Eo the characteristic adsorption 
energy of the reference material, 6 the affinity coefficient of the 
adsorbate for the solid surface, and A is the adsorption potential 


given by 

Po 
A=-—AG= RTIn (2) (11) 
Combining Eqs. (10) and (11) and then rearranging gives: 


RP? of Pi 
ln V = In Vo = BEo In P (12) 


from which a plot of In V versus In*(Po/P) should give a straight 
line, the Y-intercept of which is In Vo, which allows the limit- 
ing micro-pore volume for each sample to be determined after 
converting volume of gas adsorbed to liquid volume. Fig. 12(a) 
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Fig. 12. (a) Dubinin—Radushkevich method of micro-pore volume analysis and 
(b) changes in EMD micro-pore volume with heat treatment. 
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shows a plot of Eq. (12) for the 400 °C heat treated sample, while 
Fig. 12(b) shows the changes in limiting micro-pore volume with 
heat treatment temperature. What is apparent firstly in Fig. 12(a) 
is the only partial linearity of the complete data set. Evidently at 
both low and high In?(Po/P) values there is substantial deviation 
from the model. These regions of the data correspond to high 
and low partial pressure adsorption data, respectively, which for 
many other adsorption isotherm models, are regions were model 
breakdown also occurs. For the DR model this non-linearity has 
been noted before [19], with the significance of the results being 
subject to considerable debate. Nevertheless, using our best esti- 
mate of the linear region in the plot of Eq. (12) to generate the 
data in Fig. 12(b), we can see that the limiting micro-pore vol- 
ume follows a trend that is remarkably similar to the changes 
in BET surface area for our samples. Assuming that the empir- 
ical nature of the DR method is legitimate for these samples, 
then we can firstly confirm that the bulk of EMD surface area 
is indeed associated with micro-porosity. Secondly, the data in 
Fig. 12(b) again supports our collapsing pore model for EMD 
heat treatment. 

The Dubinin—Astakhov (DA) method [19,23] of micro-pore 
analysis is very similar to the DR method, except an additional 
degree of freedom is added to the volume filling expression; i.e. 


A n 
rower Gi) o 


Here the value of n is selected to give the best fit of the model 
to the data, with again Vo being the analysis result. As has been 
reported previously [19], the value of n was close to 2 for each 
sample, with the end result being that the micro-pore volume 
followed the same trends as those obtained for the DR method 
(Fig. 12(b)). 


3.3.4. Density functional theory (DFT) based models 

As we have seen above, conventional methods of determining 
a pore size distribution from gas adsorption and desorption data 
are quite empirical in nature. This is particularly true for micro- 
porous materials where pore sizes are comparable to adsorbate 
dimensions. A more recent advance [19] is to model the distri- 
bution of gas molecules in the system consisting of a gas and 
solid adsorbent using first principles; i.e., molecular dynamics 
or the Monte-Carlo methods. Another, less calculation inten- 
sive method involves using density functional theory (DFT) to 
describe this distribution. Despite being less demanding from 
a calculation perspective, the use of DFT provides an excellent 
representation of such an inhomogeneous system. Micromerit- 
ics DFT Plus v2.00 is an example of software that employs 
the DFT methodology. Essentially the software has a library 
of standard isotherms that have been calculated using the DFT 
methodology. These standard isotherms are specific for vari- 
ous gas-solid couples, as well as for certain pore shapes and 
sizes. The software then adjusts the relative contributions of 
these standard isotherms to fit to the experimental data, hence 
giving rise to an pore size distribution. The pore size distribution 
using such an approach, assuming slit shaped pores, is shown in 
Fig. 13. 


70 


60 


50 


40 


30 


Incremental Pore Volume / 10* cm? g' 


l 10 100 1000 10000 
Pore Size / Å 


Fig. 13. Pore size distribution of the heat treated EMD series as determined 
using a DFT based method. 


What is immediately obvious with this pore size distribu- 
tion data is the enhanced resolution that is available, covering 
around three orders of magnitude of pore sizes. In terms of 
what the data represents, it is again consistent with our pro- 
posed model for structural and morphological change. At very 
small pore sizes (~6 and ~15 A) heat treatment causes a con- 
siderable decrease in the adsorbed volume associated with these 
pore sizes. The pore volume associated with the slightly larger 
30A pores is interesting in that it remains relatively constant 
with heat treatment, except at the highest temperature (400 °C) 
where it seems as though a number of these pores have been 
closed down. At larger pore sizes (>50 A) there is a steady 
increase in pore volume due to stress-induced cracking as a 
result of material density differences, as has been discussed 
previously. 


4. Summary and conclusions 


A series of heat-treated EMD samples (100° to 400 °C) have 
been prepared and characterized using a collection of physical 
methods to evaluate the structural and morphological changes 
that have occurred as a result of heat treatment. The key conclu- 
sions that can be made from this work are: 


(i) The structure of the starting EMD (y-MnO2) undergoes 
structural contraction along both the a and b axes, but 
expansion along the c axis. This has been explained in terms 
of the thermal diffusion of manganese ions (both Mn(IV) 
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and Mn(IIJ)) resulting in a denser, less defective B-MnO2 
structure. 

(ii) As a result of this manganese ion movement under heat 
treatment, micro-pores within the solid are sintered shut, 
the net result of which is a considerable decrease in the 
surface area. 

(iii) The structural transition on heat treatment from y- to B- 
MnO, also brings about an increase in material density. 
The structural contraction that occurs as a result has been 
found to induce cracking which causes an increase in the 
number of larger pores. 
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